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Abstract. We consider the interaction between perturbations in the infiaton and in 
the metric during the preheating phase in simple inflationary models. By numerically 
integrating the Einstein field equations we are able to gauge the impact of non-linear 
gravitational effects on preheating for the first time. In the A<^ 4 model we find 
a large increase in the amplitude of sub-Hubble metric modes, beyond that due 
to gravitational collapse alone. There is significant mode-mode coupling and the 
amplification is eventually terminated by back reaction effects. We suggest that such 
enhancement of inhomogeneity will change the behaviour of the post-inflationary 
universe. 

PACS numbers: 98.80Cq, 04.25.Dm BROWN-HET-1181, TA-573 



1. Introduction 

Inflationary cosmology proposes a period of nearly exponential growth in the 
primordial universe and in so doing solves a number of the problems that face big 
bang models. In most implementations of inflation, the post-inflationary universe is 
a very cold and very empty place, and the first question one has is, where did the all 
matter seen today came from? The short answer is that it was due to a process called 
reheating in which the infiaton decayed and released its energy into "normal" matter 
and radiation. It was realised in 1990, however, that the long answer, the details of 
reheating, might turn out to be very interesting indeed. 

Traschen and Brandenbergerji] showed that if the infiaton field oscillates 
coherently at the end of inflation, then it may drive a period of explosive particle 
production through the mechanism of parametric resonance. This process is often 
termed preheating and its ramifications for cosmology include the areas of non-thermal 
phase transitions and defect production, baryogenesis, dark matter, gravitational 
waves and primordial black holes (see [|[ ||, ^] and refs. within). In models which allow 
preheating, reheating is now understood as preheating followed by thermalisation. 

In the standard analysis of preheating one assumes a Friedmann-Lemaitre- 
Robertson- Walker (FLRW) universe and then considers, both analytically^, || and 
numerically]^, |tJ, the field dynamics in this background. In general one finds that 
certain modes in the momentum space of the fields are exponentially populated, with 
the precise structure of these resonance bands being highly model-dependent. The 
result is that the fields soon deviate from homogeneity and there are significant non- 
linear effects such as rescattering of the particles, possible symmetry restoration and 
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back reaction on the inflaton. Although the FLRW scale factor can be calculated 
self-consistentlyp, |9|] in these scenarios, it was realised recently that the impact of 
inhomogeneityjig, [H], g [R§ [L4[ i HI @ and 

non-linearity || [H| [l6) on the metric 

must also be considered. 

Our analysis of preheating is new in that we use the full Einstein field equations. 
This enables us to study non-linear gravitational effects that are not accessible in 
the perturbative approximation. We find that parametric resonance can induce large 
metric inhomogeneities and that mode- mode coupling and back reaction is significant. 
At present, our principal simplifying assumption is that the inhomogeneity is in one 
spatial direction only. This means we have a f + 1-dimensional system of partial 
differential equations for two metric functions. 

The larger goal we are pursuing is to determine the spectrum of amplified 
fluctuations, in both the fields and the metric, that arise from preheating. At issue 
is whether predictions of preheating are in accord with cosmological observations. 
Recently it has been argued [|l3| [ll| that super-Hubble modes can be amplified during 
preheating. The rationale is that inflation leaves the oscillating inflaton coherent on 
scales much larger than Hubble radius. If this is the case, one might have to reconsider 
the standard perturbative account of the fluctuations that lead to large-scale structure 
formation, especially if mode-mode coupling occurs. On the other hand, it is certainly 
the case that amplified fluctuations on Hubble and sub-Hubble scales will give rise to 
gravitational waves JlT], [Tq] , and one would also like to know if primordial black hole 
formation after inflation is affected fll3|, [Tq] . 

Here we briefly survey our progress to date. We have considered preheating 
in two models of chaotic inflation, supposing there are no fields other than the 
inflaton. We have confirmed Q an analytical result |Q that the m 2 4> 2 model does 
not undergo parametric resonance. In the \<p 4 model where the standard analysis 
predicts parametric resonance in the field we have found strong amplification of 
the corresponding metric modes, and coupling of these modes to modes outside the 
resonance band. The modes which are amplified are sub-Hubble modes. We will 
report more fully on these results in a forthcoming paper|i6[. 



2. Metric and initial conditions 

For simplicity we have assumed a universe with planar symmetry [j| so the metric 
functions depend only on time and one spatial co-ordinate. We studied preheating 
after m 2 (j) 2 inflation|| using the metric 

ds 2 = dt 2 - A 2 (t, z) dz 2 - B 2 (t, z) (dx 2 + dy 2 ) . (1) 

In the limit of small spatial inhomogeneity this resembles a flat FLRW metric written 
in terms of physical time. We could also use this metric for the A</> 4 case, but it turns 
out to be more convenient to introduce a "conformal-likc" metric. This is obtained 
via a co-ordinate transformation: r\ = i](t, z) and ^ = z), which allows (|l|) to be 
written as 

ds 2 = a 2 (r?, C) (d V 2 - d( 2 ) - P 2 ( V , C) (dx 2 + dy 2 ) . (2) 

Here, we will concentrate on this second metric as we have discussed ([!]) elsewhere || 
and the way we choose the initial conditions is the same in both cases. 

f One may think of this as a universe with a flat, 2-dimcnsional, maximally symmetric subspace; in 
contrast the FLRW metric has a 3-dimensional, maximally symmetric subspace. 
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Figure 1. The Fourier transform of the metric perturbation, for a resonant mode 
is shown: the left panel gives the perturbative result, while the right panel shows the 
evolution of the mode derived from the full non-linear analysis. 



The Einstein field equations, G^ v = —kT^ v , give equations of motion, i.e. 
equations involving second derivatives in time, for a and 0. In addition one obtains 
two constraint equations which contain only first derivatives in time. The equation 
of motion for follows from T^ v '. v = 0. (We refer the reader to |l6| for the explicit 
form of these rather lengthy equations.) We numerically integrate the equations of 
motion using the same techniques of || and, in addition, simultaneously compute 
the co-ordinate transformation which takes us back to metric ([!]). This allows us to 
determine A and B as well. 

Initial conditions must be chosen to satisfy the constraint equations. This is, 
in general, quite difficult to do, but we focus on a relatively simple scenario where 
a single mode only is initially excited. Then it is possible to let 0(0, £) = 0o and 
a(0, C) = /3(0,C) = 1< Inhomogeneity is put into the system via the inflaton's 
momentum 



0,^(0,0 = 00 



2tt/cC 
Z 



(3) 



where k is the number of times the fluctuation will fit into our simulation "box" of 
length Z. The size of the perturbation is governed by e and is not required to be 
small. We start our simulations at the end of inflation and this dictates the values of 
0o a- n d 0o- Our choice of rescaling makes the initial density perturbation of order e 2 . 
The constraint equations then determine that 
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0,u(O,C) 




(4) 
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where (• ■ •) denotes a spatial average. In the limit of no perturbation these initial 
conditions ensure the metric reverts to a FLRW type. 
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Figure 2. The Fourier transform of $ for a resonant mode (larger amplitude) and 
its second harmonic (smaller amplitude), which is outside the resonance band. It 
grows because, in the non-linear analysis, it is coupled to the resonant mode. 

3. Results and discussion 

In the standard analysis, where one ignores metric perturbations, there is no 
parametric resonance in the m 2 (j> 2 model. We re-examined this model using our new 
approach and were able to show that this null result continues to hold. This was in 
agreement with a perturbative analysisjl4| which served as a check on our method. In 
addition, we began to explore the non-linear regime by making the initial perturbations 
artificially large. Now mode-mode coupling becomes noticeable and we conjectured 
that the effect would be to broaden the instability bands in a resonant system. We 
have been able to confirm this conjecture in our latest work|jl6|. 

The A</> 4 model is the simplest model which exhibits parametric resonance and 
is therefore an important test case for our approach. Tantalisingly, the modes in 
resonance are initially just inside the Hubble radius. The Hubble radius grows 
subsequently, so we only observe sub-Hubble mode amplification in this model. 

Fig. 1 compares the amplification of a resonant mode in the perturbative and in the 
non-linear analysis. To facilitate the comparison we plot the gauge invariant metric 
perturbation $[[18), which we extract from our numerical data. The first thing to 
notice is that the perturbative analysis of a resonant system always fails; the resonant 
growth takes the system away from the perturbative regime. The full analysis reveals 
what really happens: resonance terminates due to back reaction effects. An example 
of the coupling between the modes that the perturbative analysis misses can be seen 
in Fig. 2. 

Perhaps the most dramatic ramification of resonance can be see in the evolution 
of the metric component A. When we simulated the effect of an initial perturbation 
with wavelength not in the resonance band, we found only a small increase in the 
inhomogeneity due to collapse of the over-dense regions. In the left panel of Fig. 3 
one sees the result: A ~ rj, which is what one expects for a radiation dominatecjj], 
FLRW universe. The situation was quite different, however, when we considered a 
mode in the resonance band. Here the growth in inhomogeneity is initially driven by 



f The universe is initially radiation dominated in the A0 4 model. 
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Figure 3. The evolution of A (the g zz component of (|l|)) is plotted as a function of 
r) and The left panel shows the evolution of a mode with k slightly too large to be 
in resonance, while in the right panel a mode undergoes resonance, and significant 
inhomogeneity is generated. The units are arbitrary. The initial perturbation is 
e 2 = 10 — 6 , and the simulations begin at the end of inflation. 



parametric resonance and leads to a large density contrast, 8p/p ~ 1. After parametric 
resonance terminates, gravitational collapse continues to enhance the inhomogeneity. 

It is because the system is initially driven that we think the formation of primordial 
black holes might need to be reconsidered in preheating scenarios. In the usual 
understanding, one requires Sp/ p > 1 as the mode comes inside the Hubble volume, 
and black hole formation proceeds by gravitational collapse. In the At/) 4 model, even 
though the relevant metric modes are always sub-Hubble, they are being amplified 
resonantly, not gravitationally. 

We are now extending the initial results mentioned above|l6[|. First we will 
incorporate general initial conditions and then study more realistic models which 
include other fields beside the inflaton. Once we have determined the spectrum 
of amplified fluctuations we will be able to consider whether or not cosmological 
observations can be used to rule out certain preheating models. 

Acknowledgments 

We would like to thank Robert Brandenberger, Fabio Finelli, Bruce Bassett and David 
Kaiser for useful discussions. Computational work in support of this research was 
performed at the Theoretical Physics Computational Facility at Brown University. 
RE is supported by DOE contract DE-FG0291ER40688, Task A. 



References 

[1] J. Traschen and R. Brandenberger, Phys. Rev. D42, 2491 (1990). 

[2] L. Kofman, A. Linde and A. Starobinsky, Phys. Rev. D56, 3258 (1997). 

[3] M. Parry and R. Easther, Ph ys. Rev D5fl OK I 301 (1999). 

[4] A. Sornborger and M. Parry, iep-ph/981152C . 

[5] P. Greene, L. Kofman, A. Linde and A. Starobinsky, Phys. Rev. D56, 6175 (1997). 

[6] S. Khlebnikov and I. Tkachev, Phys. Rev. Lett. 77, 219 (1996). 

[7] T. Prokopec and T. Roos, Phys. Rev. D55, 3768 (1997). 



() 



[8] S. Ramsey and B. Hu, Phys. Rev. D56, 678 (1997). 

[9] D. Boyanovsky, D. Cormier, H. de Vega, R. Holman, A. Singh and M. Srednicki, Phys. Rev. 
D56, 1939 (1997). 

[10] H. Kodama and T. Hamazaki, Prog. Theor. Phys. 96, 949 (1996). 

[11] T. Hamazaki and H. Kodama, Prog. Theor. Phys. 96, 1123 (1996). 

[12] Y. Nambu and A. Taruya, Prog. Theor. Phys. 97. 83 (19 961. 

[13] B. Bassett, D. Kaiser and R. Maartens, [iep-ph/9808404 . 

[14] F. Finelli and R. Brandenberger, Phys. Rev. Lett. 82, 1362 (1999). 

[15] B. Bassett, F. Tamburini, D. Kaiser and R. Maartens, friep-ph/9901319 

[16] R. Easther and M. Parry, BROWN-HET-1152 (in preparation) (1{)99). 

[17] S. Khlebnikov and I. Tkachev, Phys. Rev. D56, 653 (1997). 

[18] V. Mukhanov, H. Feldman and R. Brandenberger, Phys. Rep. 215, 203 (1992). 



